Hysteresis and Return Point Memory in Artificial Spin Ice Systems 
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We investigate hysteresis loops and return point memory for artificial square and kagome spin 
ice systems by cycling an applied bias force and comparing microscopic effective spin configurations 
throughout the hysteresis cycle. Return point memory loss is caused by motion of individual defects 
in kagome ice or of grain boundaries in square ice. In successive cycles, return point memory is 
recovered rapidly in kagome ice. Memory is recovered more gradually in square ice due to the 
extended nature of the grain boundaries. Increasing the amount of quenched disorder increases the 
defect density but also enhances the return point memory since the defects become trapped more 
easily. 

PACS numbers: 75.10.Hk,75.60.Jk,75.10.Nr,82.70.Dd 
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Frustration effects arise in many condensed and soft 
matter systems when geometric constraints prevent col- 
lections of interacting elements such as spins or charged 
particles from simultaneously minimizing all pairwise in- 
teraction energies. One of the best known frustrated sys- 
tems are the spin ices [J, Q , named for their similarity 
to the frustrated proton ordering in water ice Spin 
ices have been realized in both two and three dimensions 
and exhibit interesting excitations such as effective mag- 
netic monopoles [J, 0] . More recently, artificial spin ices 
were created with arrays of nanomagnets colloidal 
particles [3, EBl , or vortices in nanostructured supercon- 
ductors [l7|. In artificial ices, direct visualization of the 
microscopic effective spin configurations is possible, and 
system parameters such as interaction strength, doping, 
or the amount of quenched disorder can be controlled. 
Under a varying external field, changes in the microscopic 
configurations can be imaged and used to construct hys- 
teresis loops 14 1, as shown for kagome ice where 
the motion, creation, and annihilation of topological de- 
fects along the hysteresis cycle were demonstrated 0. 
Memory effects are generally associated with hysteresis, 
and in return point memory (RPM), the system returns 
to the same microscmic configuration after completing 
a hysteresis loop l^- ll] . Recently developed techniques 
show that in real magnetic materials, RPM occurs in 
strongly disordered samples and is absent for weak dis- 
order [3, 0j El] ■ Certain classes of T = disordered 
spin systems, such as the random field Ising model [11], 
exhibit perfect RPM, while other systems require many 
loops to organize into a state with RPM 2^ 21|. 



Artificial spin ices are an ideal system for studying 
RPM since they exhibit hysteresis and the microscopic 
states can be visualized directly. The type of topological 
defect that forms and its mobility varies in different ice 
systems, ranging from mobile monopol es in kagome 
ice to less mobile grain boundaries 0, Il7l| in square ice, 
and this could modify the RPM behavior. To quan- 
tify this, we perform numerical simulations of hystere- 



sis in artificial square and kagome spin ices constructed 
from colloids in double- well traps with varied amounts of 
quenched disorder. Our model was previously shown to 
capture the behavior of square and kagome ices 1^, 1^ . 
and the number and type of topological defects present 
can be controlled by changing the amount of quenched 
disorder 17[. We using molecular dynamics simulations 
to avoid the spurious loss of RPM that occurs under 
Monte Carlo spin flip methods 22|. Our work implies 



that RPM phenomena can be studied in general artifi- 
cial spin systems and not just artificial spin ice systems, 
which would provide a new method for exploring micro- 
scopic memory effects in condensed matter systems. 

Simulation- We simulate an artificial spin ice of N 
charged colloidal particles trapped in an array of elon- 
gated double-well pinning sites that have two states 
determined by which well is occupied by the colloid. 
The dynamics of colloid i is governed by the over- 
damped equation of motion rj{d'R.i / dt) = F^'^-I-Ff -|-Fext, 
where the damping constant -q — 1. The colloid-colloid 
interaction force has a Yukawa or screened Coulomb 
form, F'^'^ = 



(l/i?y) exp(-K_Rjj )f y. Here R, 
(Ri — 'R.j)/R. 



Ri — R-i , Ri 



is the position of particle 
Fq — Z*''^ /{Aneeo), Z* is the unit of charge, e is the sol- 
vent dielectric constant, q is the dimensionless colloid 
charge, 1/k = 4ao is the screening length, and oq is the 
unit of distance that is typically of order a micron. We 
neglect hydrodynamic interactions between colloids since 
we work in the low volume fraction limit and the colloids 
remain confined in the pins. The pinning force F^ arises 
from Np elongated traps of length I — 1.333ao, width 
dp = 0.4ao, and depth fp = lOOi^o- The pin ends are 
parabolic confining potentials with radius rp = 0.2ao. A 
cylindrical force restricts motion in the direction perpen- 
dicular to the long axis of the pinning site, and a barrier 
in the center of the pinning site is produced by a repul- 
sive parabolic force of height fr that creates two energy 
minima on either end of the pin For square ice the 




FIG. 1: Schematic of a portion of the square ice sample. A 
charged colloid (dots) can sit in either end of each trap (rect- 
angles), (a) A biased state for foxt applied ai 9 = 45° from 
the X axis, where all the vertices are in the (0011) configura- 
tion, (b) The negative biased state after the drive is reversed 
contains only (1100) vertices. By cycling Fcxt between these 
two extremes we construct an effective magnetization curve. 
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pins are arranged with v = A traps meeting at each ver- 
tex, as shown in Fig. 1, while for kagome ice, u = 3 traps 
meet at each vertex 17|. The distance between adjacent 
vertices is a = 2ao and there are Ny vertices. Our square 
ice has 35 x 35 vertices (iY„ = 1225) and Np = 2450 
elongated pins, while our kagome ice has 40 x 40 vertices 
{Ny — 1600) and Np — 2400 elongated pins. Disorder is 
added to the system by increasing or decreasing fr in in- 
dividual pinning sites according to a normal distribution 
with mean fb = l.Oi^o a-nd standard deviation a. This 
is analogous to varied island coercive fields in the nano- 
magnetic system. We initialize the system by placing a 
colloid in one randomly selected end of each pinning site 
so that N = Np. To construct a hysteresis loop we apply 
an external force Fext uniformly to the sample, which for 
charged colloidal particles could be achieved using an ex- 
ternal electric field. In the kagome ice Fcxt = ^cxtx while 
in the square ice Foxt — ^cxt\/2(x -I- y)/2, as shown in 
Fig. 1. We sweep Fcxt from zero to a positive maximum 
value i^max, then back down through zero to a negative 
maximum value — J^max, and finally back up to zero to 
create one loop. 

The vertices formed by the meeting points of the pins 
are categorized by ni„ = Yli=i ^ij where c; = 1 if a col- 
loid is sitting in the end of the pin closest to the vertex 
and Ci = if the colloid is sitting in the end of the pin 
furthest from the vertex (isl . [l7] . Each possible vertex 
type can be written as (C1C2C3C4) for the square ice and 
as (C1C2C3) for the kagome ice, starting with the pinning 
site lying along the positive x axis with respect to the 
vertex and proceeding counterclockwise around the ver- 
tex. The Tiin = 2 vertices in square ice do not all have 
the same energy. We denote the ground state vertices 
(1010) and (0101) as type 2gs, and the higher energy bi- 
ased vertices (1100), (0110), (0011), and (1001) as type 
26. 

We first show that our model captures the hysteretic 
behavior observed in artificial ice systems 
the absence of quenched disorder or drive, we find ice- 



FIG. 2: (a,c,e) Square ice sample with a — 0.1. (b,d,f) 
Kagome ice sample with a = 0.1. (a,b) The reduced mag- 
netization m vs fcxt. Saturation occurs at m = ±1.0 when 
all the vertices are in biased states. Dark line: Saturated 
loop with Fmax = 2.0. Light lines: Consecutive loops with 
fmax = 0.7, below saturation. The virgin curves are not 
shown. (c,d) Fraction of unbiased vertices Nub vs _Foxt for 
the saturated loop with -Fmax ~ 2.0. (e,f) N^t vs Fcxt for re- 
peated unsaturated loops with Fmax = 0.7, with cycle number 
n increasing from top to bottom. For clarity, we omit the hor- 
izontal lines connecting Fcxt = iFmax to Fext = 0. There is a 
much greater decrease in Nub for the square ice than for the 
kagome ice. Inset of (f): min(g), the effective spin overlap in 
the n = 2 cycle vs a for (circles) kagome and (squares) square 
ice. Samples with stronger disorder have higher q values. 



rule obeying states that are ordered ground states in the 
square ice [l5| and disordered in the kagome ice. When 
we add quenched disorder with tr > 0, due to its lack of 
extensive degeneracy the square ice forms grain bound- 
aries composed of non-ice rule obeying vertices as shown 
in simulation 17 1 and experiment [6|, while in kagome 
ice isolated non-ice rule defects appear 0, Under 
the biasing drive Fext shown schematically in Fig. 1(a), 
all the square ice vertices adopt the (0011) configuration. 
For kagome ice biased along the x direction, all the ver- 
tices enter the (Oil) state. We define the reduced mag- 
netization m as the fraction of vertices aligned with the 
saturation direction, M — N.^-^ where Si = 1 for 

(0011) vertices in square ice or (Oil) vertices in kagome 
ice. Si — —1 for (1100) vertices in square ice or (110) 
vertices in kagome ice, and Si — for all other vertices. 
In Fig. 2(a,b) we plot the hysteresis loops for square and 
kagome ice samples with a = 0.1. The thick curve is 
obtained with Fmax = 2.0, beyond the saturation level 
where to = ±1. We^lot the fraction of unbiased ver- 
tices, Nub = N-^J2i=lHsi), vs fext with Fmax = 2.0 

in Fig. 2(c,d). The completely ordered biased states are 
only destroyed for 0.4 < \Fext\ < 1-0, close to the co- 
ercive fields at which the effective spin direction fiips. 
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FIG. 3: Effective spin overlap q vs F^xt during consecutive 
hysteresis loops for (a) a square ice sample with a — 0.1 and 
(b) a kagome ice sample with a = 0.1. Heavy line: Saturated 
loops with -Fmax = 2.0, including the virgin curve. Thin lines: 
Unsaturated loops with -Fmax = 0.7, with n increasing from 
bottom to top. Solid lines: clockwise loops; dashed lines: 
counterclockwise loops. In the kagome ice, q approaches 1 
after only a few cycles, while a much larger number of cycles 
are required before q approaches 1 in the square ice. 



The shape of the hysteresis loop and the peaks in the 
non-biased defect density in Figs. 2(b,d) are in excellent 
agreement with the digitally constructed hysteresis loops 
produced in experiments on nanomagnetic kagome ice 
samples 0-01 • The thin lines in Figs. 2(a,b) show con- 
secutive hysteresis loops obtained below saturation with 
^max = 0.7, near the middle of the range of i^oxt in which 
the largest number of defects appear. In Fig. 2(e,f) we 
plot Nub versus Foxt for the unsaturated hysteresis loops. 
For the square ice. Fig. 2(e) shows that N^b decreases 
with increasing n, where n is the number of loops per- 
formed, indicating that defect annihilation is occurring. 
For continued cycling beyond the number of loops shown 
in the figure, the system settles into a steady state. In 
the kagome ice. Fig. 2(f) shows that N^b hardly changes 
from one cycle to the next, indicating that only a small 
number of defects annihilate. For the saturated case with 
Fmax = 2.0 shown in Fig. 2(b,c), the Nub curves do not 
evolve under repeated looping since the sample loses all 
memory of the microscopic configuration near the coer- 
cive field once saturation is reached. 

We quantify the RPM by measuring the overlap q in 
the effective spin configurations along a hysteresis loop 
2^ 22 1 at equal values of Foxt after n complete cycles. 
For each trap, we define an effective spin S'i = 1 if the 
colloid is sitting in the right or top end of the trap, and 
Si = —1 a the colloid is sitting in the left or bottom end 
of the trap. Writing the value of Si after n cycles as S,^ 



g(Fe.0=iV-1^5l"-^)(Fe.t)5f^(Fe.t). (1) 



The term in the sum is 1 if the trap was biased in the 
same direction both before and after the complete cy- 
cle, and -1 if the colloid jumped to the other end of the 
trap. In Fig. 3(a,b) we plot q versus i^ext for both the 
saturated and unsaturated hysteresis curves in the square 
and kagome ices shown in Fig. 2(a,b). In the case of the 
saturated loops, q for the virgin curve in Fig. 3(a) shows 
that since the sample was not initialized in a biased state, 
the initial configuration differs significantly from the ef- 
fective spin configuration obtained one cycle later, but 
for n = 2 and above, q = 1, indicating perfect memory. 
For the unsaturated loops obtained with i^max = 0.7, q is 
low during the first cycle, but as n increases q gradually 
converges to a value just below q = 1. A comparison with 
Fig. 2(e) indicates that the increase in memory with in- 
creasing n is correlated with a decrease in Nub, although 
for this value of F^ax there are always some defected 
vertices present even after the system reaches a steady 
state in which the grain boundaries cease to evolve. The 
kagome ice in Fig. 3(b) shows a similar behavior except 
that q approaches 1 after only a few cycles, leading to 
a much faster establishment of RPM than in the square 
ice. In Fig. 2(f) we show that the number of defected 
vertices remains nearly constant in the kagome ice even 
under repeated cycling. This indicates that although the 
kagome ice defects do not annihilate, they are mobile 
during the first few cycles and then become pinned. Our 
results demonstrate that for the square ice, changes in 
the amount of RPM are primarily associated with the 
annihilation of defects, while in the kagome ice, RPM is 
suppressed by the motion of defects. 

Although the number of defects Nub in both types of 
ice increases with increasing disorder a, the amount of 
RPM increases with increasing disorder. We illustrate 
this in the inset of Fig. 2(f) where we plot the value of q 
on the n = 2 plateau versus a. A similar effect was ob- 
served for real magnetic systems and in spin simulations 
liL 2oL 2^ . In our system, q increases with increasing 
disorder due to the stronger pinning of the domain walls 
in the square ice or of the individual defects in the kagome 
ice. In the square ice, the disorder prevents the domain 
walls from coarsening. It was previously shown that as 
the particle-particle interaction strength in our system 
is reduced, non-ice-rule obeying vertices begin to appear 
since their energetic cost decreases [3] ■ For noninteract- 
ing colloids, the sample is strongly disordered but also has 
perfect RPM since the defected configurations arc con- 
trolled only by the local disorder and are not modified by 
particle interactions. Thus we expect that in the experi- 
mental nanomagnetic artificial ices, when the coupling is 
reduced for increased spacing between the nanomagnets. 
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FIG. 4: Vertices in (a) square ice and (b) kagome ice samples 
from Fig. 2 during repeated hysteresis cycles. Vertices are 
colored depending on how many cycles the vertex spends in 
a defect state, ranging from dark blue for never defected sites 
to dark red for permanently defected. In (a), motion and 
annihilation of defects occur near grain boundaries. In (b), 
individual defects move and are pinned independently without 
forming grain boundaries. 



the system should show increased or perfect RPM. 

To illustrate the defect dynamics, in Fig. 4 we plot the 
vertices colored according to the number of hysteresis 
cycles each vertex spent as a defected site. Red vertices 
indicate locations where defects became trapped. Fig- 
ure 4(a) shows that in the square ice, the defects orga- 
nize into grain boundaries which move and coarsen under 
repeated hysteresis cycles. In Fig. 4(b), the kagome ice 
contains no grain boundaries but has a smaller fraction 
of intermediately colored vertices compared to the square 
ice since the isolated defects become trapped after only 
a few cycles. The square ice grain boundaries are less 
well pinned than the isolated kagome defects since they 
are extended objects. The motion of individual defects 
in kagome ice has already been imaged in experiments; 
it would be interesting to observe whether these defects 
become localized within a few hysteresis cycles. 

For most artificial ice systems, thermal effects are not 
relevant; however, thermal fiuctuations can be significant 
in a colloidal system. We find that our results are robust 
against the addition of weak thermal disorder, and that 
for T > there is only a slight reduction in the asymp- 
totic value of q and a slight increase in the number of 
cycles required to reach a steady state. For higher tem- 
peratures, RPM is lost even when the system is cycled to 
saturation since the thermal fiuctuations cause random 
effective spin flips that change the path on each cycle. 
There is also no increase in RPM under repeated cycles 
at higher temperature [23|. 

In summary, we have studied hysteresis and return 
point memory effects for artificial square and kagome ices 
at the microscopic level. In the square ice for repeated 
unsaturated hysteresis loop cycles that extend to biases 
near the coercive field, the RPM increases with each cy- 
cle as the grain boundaries present in the sample coarsen 



and become pinned. In kagome ice the number of defects 
remains nearly constant under repeated hysteresis cycles 
and there is much higher RPM. Here, individual defects 
hop rather than annihilating and are eventually pinned 
at sites with stronger disorder. The grain boundaries in 
the square ice are more mobile than the individual defects 
in the kagome ice since they are extended objects. Our 
results can be tested readily in different types of artificial 
ices and also could be studied in more general artificial 
spin systems. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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